Tradeoff between Leakage and Dephasing Errors in the Fluxonium Qubit 
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We present a tradeoff between leakage and pure dephasing errors for the fluxonium circuit. We 
show that in the insulating regime, i.e., when the persistent current flowing across the circuit is 
suppressed, the pure dephasing rate decreases exponentially as the impedance of the circuit is 
increased. In contrast to this exponential decrease, the qubit remains sufficiently anharmonic so 
that gates times can still be short, allowing for significant reduction in the computational error rates. 



Superconducting circuits are a promising technology 
for building scalable quantum processors in the solid- 
state pp. They feature a tunable anharmonic spectrum 
that can be used as a two- level system, defining a qubit. 
Single- and two-qubit gates have been demonstrated [2]- 
3] and the lithographic fabrication process is a mature 
technique, allowing for large scale production. However, 
as a general rule, solid state qubits couple much more 
strongly to their surrounding environments than systems 
such as atoms or photons, leading to far shorter decoher- 
ence times. Considerable effort has thus been invested 
in fighting the causes of decoherence in different super- 
conducting circuits [SHU], in order to attempt to bring 
achievable error rates close to fault- tolerant thresholds. 
Importantly, recent experimental advances point to pure 
dephasing (i.e. dephasing without population relaxation) 
as a significant source of decoherence in superconduct- 
ing qubits coupled to a cavity [T2HT4] . In particular, in 
cases where relaxation processes can be greatly reduced 
[10J [14j [19] , pure dephasing can lead to decoherence (T2) 
times that are significantly shorter than the lifetime limit 
of2Ti. 

Of course, this need not be a problem if quantum gates 
can be performed sufficiently quickly, so that the error 
rate per gate is low. Given a dephasing rate T^, and in 
the absence of relaxation processes, the dephasing error 
probability per gate is p^ ex 1^ x r, where r ex 
is the time taken to perform a gate, and Q is the Rabi 
frequency associated to the driving field. It is thus crucial 
for low error-rate operation that the ratio of to Q can 
be made small. One way to achieve this would be to 
drive the system strongly, so that gates are performed 
extremely quickly on the timescale set by 1^. However, 
this is ultimately limited by the anharmonicity 5 offered 
by the system, i.e., the energy difference between the first 
and the second transition energies, 5 — (A21 — A10). A 
driving field resonant with the transition |0) |1) will 
induce leakage into level |2) with probability 
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This probability has to remain well below the threshold 



for loss errors [15j [16] to guarantee an error-free compu- 
tation at the logical level. Hence, a strategy for suppress- 
ing the system pure dephasing rate without significantly 
compromising its anharmonicity is attractive, as it would 
allow the error probabilities due to both pure dephasing 
and leakage to be made small simultaneously. 

Here, we show that by appropriate choice of operat- 
ing regime it is possible to exponentially increase pure 
dephasing times in a superconducting circuit, while re- 
ducing its anharmonicity only at a smaller (algebraic) 
rate. We thus obtain a tradeoff between the anhar- 
monicity, related to leakage probability, and computa- 
tional errors caused by pure dephasing. In Ref. [7] a 
similar issue was considered for the transmon circuit, 
where it was found that operating the charge qubit at 
large values of the Josephson energy provides protec- 
tion against charge fluctuations, without compromising 
the qubit controllability. In contrast, we address this 
point for the fluxonium circuit [T7] , which is a partic- 
ularly interesting superconducting circuit exhibiting ro- 
bustness against low frequency variations of both offset 
charge and external flux. Fig. [IJa) shows a simplified 
version of the fluxonium circuit, which is operated in the 
regime E L /E C <C Ej/E c ~ 1, where E L , E c , and Ej 
are the inductive, capacitive, and Josephson energies, re- 
spectively. 

Closing a superconducting island with an inductor 
breaks the confinement of charge carriers, and as a con- 
sequence the conjugate charge in the island is not neces- 
sarily an integer number of Cooper pairs that tunneled 
in or out of the island. Equivalently, since periodicity 
of the superconducting phase implies that the conjugate 
charge must be discrete, in this case the superconducting 
phase <j) is not periodic. Protection against low frequency 
charge noise is obtained for the fluxonium qubit in the 
same way as in other persistent current circuits, namely, 
by shunting the island with an inductive element in such 
a way that offset charge can be accounted for by choos- 
ing a different gauge [17]. Also, whenever the persistent 
current flowing across the circuit is suppressed, its depen- 
dence on external flux fluctuations is also highly reduced. 

Together with an enhanced protection against low fre- 
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FIG. 1. (a) Fluxonium circuit diagram. The regime at which 
the device is operated corresponds to El/Ec <C Ej/Ec ~ 1. 
The gate charge n g is not depicted since the energy spec- 
trum does not depend on it, but rather on its temporal vari- 
ation, (b) Potential V{<j)) — El4> 2 — Ejcoscj) (dashed) and 
ground state wavefunction (solid) for the values Ej/Ec = 1 
and E L /Ec = 10" 4 . 



quency noise, fluxonium also features a suppression of 
quantum fluctuations of charge. As we shall shortly see, 
quantum fluctuations of flux are proportional to the cir- 
cuit's impedance Zq = V Ec ^ lR Q_^ w j iere — h/(2e) 2 
is the resistance quantum. In the ground state, charge 
fluctuations will be inversely proportional to flux fluctu- 
ations, so in order to keep fluctuations well below one 
Cooper pair, the impedance must be made much larger 
than the resistance quantum Rq [18] . In order to attain 
this, elements with large L need to be created without in- 
creasing the associated capacitive energy of the inductive 
element, since the phase-slip amplitude decreases expo- 
nentially in the ratio y/Ej/Ec [see Eq. Q]. This has 
already been achieved in some experiments using chains 
of Josephson junctions [181120]. Of special interest is 
Ref. [20 , where the attainable ratio \JEq /El is pre- 
dicted to surpass the resistance quantum by more than 
two orders of magnitude. Moreover, recent experiments 
show that embedding a circuit within a resonator or a 
3D cavity increases its relaxation times due to screening 
of the vacuum impedance, making pure dephasing a sig- 
nificant source of decoherence [TO] [14] [19] . Also, dephas- 
ing due to measurement back-action cannot be overcome 
completely [21]. Our result is relevant for these scenarios. 

Decoherence due to coupling of the qubit to the in- 



ternal degrees of freedom of the superconducting ma- 
terial still poses a problem [22] [23]. A linear (ohmic) 
dissipative element can be modelled as an infinite set of 
non-interacting harmonic oscillators coupled to the cir- 
cuit 24 . The total Hamiltonian is 



H = E c n 2 + E L ^ 



- Ej cos( 
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2ct>Y J E L ,<Pu (2) 



where Eq 
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2C 



is the capacitive energy of the circuit, 



El = is the energy associated with the linear 

inductive element, and Ej = is the Josephson en- 
ergy. The superconducting phase across the Josephson 
junction, 0, and the number of Cooper pairs stored in 
the island, n, are quantised such that [0, n] = z, while 

Q _ 2_jtO - g ^ e dimensionless flux threading through 

i — i 
the circuit. The first line in Eq. ([21) corresponds to the 

fluxonium Hamiltonian, whose numerical diagonalization 
without bath coupling is depicted in Fig. [TJb). The sec- 
ond line gives the environment and interaction Hamil- 
tonians (subscripts i refer to environmental degrees of 
freedom), from which a master equation is derived. 

Varying the ratio Ej / Eq-, using for instance a SQUID- 
like Josephson junction, dramatically changes the spec- 
trum of fluxonium, and allows insight into the conditions 
under which a suppression of pure dephasing is expected. 
For small values of Ej/Ec, phase slips are common and 
the superconducting phase does not have a well-defined 
value. The current flowing through the circuit is negli- 
gible - it is in an insulating phase. In this regime, the 
lower-lying eigenstates of the circuit can be well approx- 
imated by harmonic oscillator states with an effective 
transition frequency uo* = v | c * , where Ec* = 2C^' 

The effective capacitance C* is given by 



C* = (2e) 2 



d 2 e(n 
dn 2 



(3) 



n=0 



where e(h) is the charge dispersion relation, and 2eh 
is the charge analogy of the quasimomentum in a one- 
dimensional lattice. For Ej/Ec > 1, the tight-binding 
approximation applies and the charge dispersion relation 
can be written as e(n) = — 2tcos(27rn) [see Fig. [2|a)]. 
The tunneling amplitude t can be obtained via instan- 
ton methods [25] , which with our energy definitions is 
given by the expression 



t 



4(2E 3 jE c ) 1/4 



exp 



Ej 
2E C 



(4) 



It is thus possible to derive an analytic expression for the 
effective capacitance, using Eq. Q and the tight-binding 
approximation, to give 



C* 



(2e) 2 
8tir 2 



(5) 
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FIG. 2. (a) Dispersion relation e(n) of the circuit for the 
ground state (blue) and the first excited state (red) for 
Ej/Ec — 0.5. Close to the h — the ground state energy 
is almost parabolic with a modified curvature, caused by the 
presence of an anticrossing at h — ±1. (b) Effective capaci- 
tance C* / (2e) 2 , in units of the capacitive energy Ec- The red 
(solid) curve corresponds to C* computed numerically using 
Eq. ([3b, while the dash-dotted line is calculated from Eqs. Q 
and (pi. (c) Variance of the wavefunction at — for a 
range of values, Ej/E c G {0, 1/5, 2/5, 2}, where Ej/E c 
increases in steps of 1/5 from top to bottom. The solid lines 
are obtained by numerically diagonalising the fluxonium cir- 
cuit Hamiltonian, and the circles correspond to the variance 
derived using the effective capa citance p icture. Note that the 
approximation gets better as yj Ec/El increases. 



which can in turn be used to obtain the variance of the 
wavefunctions through the usual harmonic oscillator re- 
lation a\ = ^^ff • Figs. j^b) and (c) depict the 
effective capacitance of the circuit computed both nu- 
merically and analytically, and the fluctuations of the 
phase as a function of the ratio ^/Ec/El, respectively. 

The effective capacitance description will break down, 
however, when the wavefunction in the superconducting 
phase basis does not expand over several minima, i.e., 
the phase particle does not "see" a crystal anymore. In 
the charge basis, this happens when the wavefunction is 
so wide that that the curvature of the dispersion relation 
at the origin h = 0, which arises in the expression for the 
effective capacitance, cannot be evaluated. We expect 
this to happen when the wavefunction starts to become 
localized in the phase basis, so a persistent current can 
be measured. 
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FIG. 3. Maximum value of the persistent current, normalised 
by the inductive energy, I p /El- Increasing the ratio Ej/Ec 
causes the circuit to enter a superconducting regime. This 
transition depends also on the ratio El/Ec, since smaller 
inductive energies result in a wider wavefunction, and the 
transition then takes place for larger values of the Josephson 
energy. The figures on the left illustrate the shape of the 
wavefunction at different points of the diagram. Whereas the 
phase is distributed roughly according to a Gaussian distri- 
bution in the insulating regime, it becomes localised in the 
superconducting regime. 



As Ej/Eq increases, phase slips are exponentially re- 
duced [Eq. Q] and the phase does indeed become lo- 
calized. In this regime a persistent current arises and 
the ground state strongly depends on the external flux 
6. The transition from the insulating to the supercon- 
ducting regime is illustrated in Fig. [3] In the insulat- 
ing regime, the wavefunctions spread over several min- 
ima of the periodic potential and fluctuations of the ex- 
ternal flux, which cause the relative offset between the 
parabolic and periodic components of the potential to 
fluctuate, hardly influence the wavefunctions. In the su- 
perconducting regime, the wavefunctions are localized in 
a minimum of the periodic potential and are highly sus- 
ceptible to fluctuations of the external flux. The results 
of Fig. [3J and the preceding analysis, therefore suggests 
that in order to effectively suppress pure dephasing we 
should look at regimes of large Ec/El, as this will allow 
us to remain in the insulating regime for a larger range of 
Ej/Ec, thus minimising the adverse effect on the circuit 
anharmonicity. 

Writing a master equation for the circuit truncated to 
its two lowest-lying eigenstates allows us to identify the 
pure dephasing rate of the qubit, oc = |(1|0|1) — 
(0|(/>|0)| 2 . Our numerical simulations show that in the 
insulating regime pure dephasing is in fact exponentially 
suppressed as ^ Ec/El increases, as shown in Fig.^a). 
This is consistent with protection against dephasing for 
low frequency noise [17], since variations of the energy 
difference between the two lowest lying eigenstates are 
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FIG. 4. (a) The dephasing rates are proportional to = 
|(1 10| 1) - (O|0|O)| 2 . The different lines correspond to the 
values Ej/Ec G {1/5, 2/5, 1}, where Ej/Ec increases in 
steps of 1/5 from bottom to top. In the insulating regime, the 
suppression can be of more than ten orders of magnitude, (b) 
With the same colour code as above (Ej/Ec now increases 
in steps of 1/5 from top to bottom), this plot shows that the 
relative anharmonicity 5 r decays only algebraically. This an- 
harmonicity will ultimately determine the minimum length of 
the pulses. 



exponentially suppressed in ctq = \J Ec* /El: 
^A 10 = -Ej± «1| cosW> - 6>)|1> - <0| cos(0 - 0)|O» , 

*-Ej8m(0)(T%e-$, (6) 

assuming that the effective capacitance approximation 
holds so that the wavefunctions are approximately har- 
monic. 

Turning now to the effect on anharmonicity, a useful 
parameter to characterise the quality of the two-level sys- 
tem truncation is the relative anharmonicity, 5 r = 5/Aiq. 
As shown in Fig. Qb), our calculations predict that the 
relative anharmonicity decays only algebraically as the 
impedance increases in the considered regime, which is 
in stark contrast with the exponential reduction of the 
dephasing matrix element shown in Fig. Qa). As can be 
inferred from Eq. ([!]) , in order to keep the leakage proba- 
bility small, pl <C 1, the gate time r must remain longer 
than the inverse anharmonicity, i.e., r > h/S. Shortest 
pulse lengths currently achievable are of the order of 10 
ns [10], which remains longer than K/ (5 r Aio) in the range 
of values of S r we obtained, and for typical transition fre- 
quencies of about 2tt x 10 GHz [T9] . 

Finally, in Fig. [5j we show a tradeoff between the rel- 
ative anharmonicity and the pure dephasing matrix el- 
ement M^, related to leakage and dephasing errors, re- 
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FIG. 5. Tradeoff between pure dephasing matrix element 
and the relative anharmonicity, obtained from combining the 
plots in Figs. [4^a) and (b). The various curves correspond 
to the values Ej /Ec £ {1/5, 2/5, 1}, increasing in steps of 
1/5 from top to bottom. For Ej/Ec = 2/5 and above, the 
relative anharmonicity is large enough to ensure low leakage 
rates even at very high impedances. 

spectively. We observe that, in the insulating regime, it is 
possible to exponentially reduce the dephasing probabil- 
ity per gate without compromising significantly the tol- 
erance to leakage errors. The threshold for loss-tolerant 
quantum computing is significantly higher than for com- 
putational errors [16 , and leakage out of the computa- 
tional subspace is an inconvenience relatively easier to 
circumvent than detecting and correcting unknown com- 
putational errors. Whenever possible, it is thus desirable 
to trade computational errors for loss errors, which can 
be done by increasing the impedance of the circuit. 

Summary - We have shown that it is possible to expo- 
nentially suppress the pure dephasing rate of the fluxo- 
nium qubit caused by the presence of a Markovian bath, 
at the expense of (slightly) reduced anharmonicity. The 
anharmonicity only decays algebraically, however, allow- 
ing for a tradeoff in which the length of the driving pulses 
need not be increased significantly in order to compen- 
sate for leakage errors. Thus, both the dephasing and 
leakage error rates may be made small simultaneously. 
At larger values of the circuit's impedance, the lowest- 
lying states are better approximated by harmonic oscil- 
lator states in the insulating regime, and the transition 
from the insulating to the superconducting phases occurs 
at increasing values of the ratio Ej/Eq- The indepen- 
dence of the harmonic oscillator states to external flux 
variations can thus be inherited by the fluxonium states, 
though without a complete loss of anharmonicity nec- 
essary to define a qubit. Interestingly, this is achieved 
in the same regime in which low-frequency external flux 
noise is suppressed. Experimental advances towards re- 
alisation of very high impedances [20] , as well as towards 
to the inhibition of relaxation processes in 3D cavities 
[T0l H4] . make our result relevant in the design of the 
next generation of superconducting qubit s. 

This work was completed shortly after the tragic death 
of one of the authors, Sean Barrett. It is with immense 
gratitude that we acknowledge his close friendship, his 
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